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Abstract 



We studied the geometrical and topological rules underlying the dispositions and the size distribution of non- 
overlapping, polydisperse circle-packings. We found that the size distribution of circles that densely cover a plane 
follows the power law: N(R) oc R~ a . We obtained an approximate expression which relates the exponent a to the 
average coordination number and to the packing strategy. In the case of disordered packings (where the circles 
have random sizes and positions) we found the upper bound a Max = 2. The results obtained for circles-packing 
was extended to packing of spheres and hyper-spheres in spaces of arbitrary dimension D. We found that the size 
distribution of dense packed polydisperse D-spheres, follows -as in the two dimensional case- a power law, where 
the exponent a depends on the packing strategy. In particular, in the case of disordered packing, we obtained 
the upper bound uuax = D. Circle-covering generated by computer simulations, gives size distributions that 
are in agreement with these analytical predictions. Tin drops generated by vapour deposition on a hot substrate 
form breath figures where the drop-size distributions are power laws with exponent a ~ 2. We pointed out the 
similarity between these structures and the circle-packings. Despite the complicated mechanism of formation of 
these structures, we showed that it is possible to describe the drops arrangements, the size distribution and the 
evolution at constant coverage, in term of maximum packing of circles regulated by coalescence. 



1 Introduction 

The plane is not very efficiently covered with circles, but circle-covering is a good model-system for the study of 
the formation and evolution of many natural and artificial systems as granular materials 111 |^] , island formation in 
metal films , segregation problems || , plate tectonics and turbulence || [7) , focal arrangement in smetic liquid 
crystals |fl8|, |l9|. One interesting problem related with packing circles is to find the size distribution that efficiently 
leads to the densest covering, compatible with the packing strategy. In some particular cases, as the Apollonian 
covering j2Cj (where the circles are packed all tangent to each other following a defined sequence), analytical and 
numerical solution are available in literature ^, ^, ^). In this paper the problem is discuss in the general case 
of disordered packing (or "osculatory" packing), where the circles are set with random sizes and positions. The 
motivation of this work is understand the basic mechanisms which are underlying the morphogenesys of breath 
figures [jl0| |ll[ |l2|, |l^, |lj, [l6| . Water which condense on a glass form a densely-packed system of droplets known 
as breath figures. The formation of the droplet system is regulated by two basic mechanisms: the independent 
grow of a drop supported by the condensation and the melting of two or more drops in a coalescence phenomena. 
These two mechanisms leads to systems with a wide range in the drop sizes. The distribution is characterized by 
rather uniform drop sizes in the region of large radii (drops which are grown from the originally firstly nucleated 
drops trough a chain of independent grow and coalescence) and a power law in the region of small radii (drops 
which are generated by re-nucleation on the surface liberated by coalescence). The a formation of these structures 
is an extremely complicate dynamical process. The point of view here adopted is that the system-morphogenesys 
is mainly driven by the geometrical constrains which regulate the drop-packing. 

Two dimensional structures generated by packing circles (or by packing others isotropic natural objects like dew 



drops on a glass), have some properties which are independent of the specific formation mechanism. The main 
similarities can be schematized in the following three points: (i) wide range in the sizes; (ii) scale invariance in 
the packing arrangement; (iii) power law in the size-distribution (N(R) oc R~ a ). Let us briefly discuss about the 
possible physical origin of these three similarities: 

i) A wide range in the particle sizes is a necessary condition to reach dense packing. For example, the packing 
of circles with equal radii Rq has a maximum density of 0.907. The density can be raised to p — 0.95 by filling 
the interstitial spaces with circles of sizes R\ = Rq/6A. Another increment to p — 0.97, requires interstitial 
radii equal to Ri — Rq/15.9. In general, the density can be increased up to any value p < 1, but this requires 
the utilization of interstitial circles with dimensions which rapidly decrease. 

ii) Consider a procedure where a dense circle-packing is generated by filling the interstitial spaces with circles 
of maximum sizes compatibly with the condition of non overlapping. (Note that, fixed the maximum radius, 
this procedure generates structures which minimizes the surface extension at fixed covered area.) In this 
procedure the only relevant metric parameter is the ratio between the external radii of the circles that 
generate an interstice and the internal radius of the circle that fill this interstice. For example, in the 
Apollonian case, this ratio tends to the value x ~ 2.9. In general, x tends to a constant value that is related 
to the packing strategy. It is straightforward to see that packing characterized by constant values of x are 
scale-invariant . 

iii) The number of circles introduced at any covering step increases following a geometrical progression. For 
example, in the Apollonian packing, if one starts with four circles in contact, at the first covering step one 
has 3 interstitial spaces to fill by circles, at the second step the interstices to fill are 9, at the third step this 
number is 27 and, at the v th are 3". In the general case, the number of circles introduced at any stage grows 
as a geometrical progression a v . By associating this increment in the number of interstitial circles with the 
condition of scale invariance, one obtains a size distribution that follows the law N(R) oc R~ a , where the 
coefficient a depends on the packing strategy. 

In this paper we obtain an approximate expression that relates the exponent a to the local topological properties 
of the circle-packing and to the packing strategy. The limiting value of such coefficient in the case of disordered 
maximum packing have been found equal to 2. 

Simple computer simulations of disordered circle covering have been performed. The results, presented in section 
^, confirm the theoretical predictions: the circle size distribution follows a power law with exponent a ~ 2 



The size distribution of tin droplets condensed on a hot, flat surface have been experimentally studied |2l|, 22 
|23| . p4, 25 1 . In section^, the experimental data are interpreted in terms of maximum circles packing. Through this 



point of view it is possible to give a simple explanation of some statistical and dynamical characteristic of these 
systems. In particular, the drop size distribution, the coverage and the self similar evolution at constant coverage, 
have been interpreted in the framework of maximum packing of drops regulated by coalescence. 



2 Topological and geometrical rules in packing 

Consider a plane densely covered by circles placed at random following the only constraint of non overlapping. 
Consider the Dodds network Q constituted by the edges that connect the centres of circles in contact (see. fig-(0) 
). Such a network has a number of vertices (in the following indicated by V) equal to the number of circles. 
The vertex connectivity ((z)) is the average number of neighbours in contact with a circle. The other network 
parameters (in particular the number of edges E and the number of faces F) are given, in term of V and (z), by 
the following relation [^7j 

(z)V = 2E (1) 
(any vertex is surrounded by (z) edges and any edge is bounded by 2 vertices^) , and by the Euler's formula 

V-E + F = X , (2) 
where \ is the Euler Poincarre characteristic and takes the value of 1 for an Euclidean plane. 

1 This relation holds only if in the network there aren't insulated vertices. This condition is generally satisfied in dense packed 
systems. 
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Figure 1: The Dodds network generated by conneting the centers of circles mutually in contact. 



The combination of equations (|lj) and (|2|) leads to a useful relation between the number of circles V and the 
number of interstitial spaces between the circles (that is equal to the number of faces (F) of the Dodds network), 
we have: 

F=(&-l)V + x ■ (3) 

Note that at the point (z) = 2 the interstitial spaces between circles become connected (F = x = 1, in the 
Euclidean plane) and, consequently, the network unconnected. This is the percolation threshold in two dimension. 
For the purposes of the present paper it is interesting the opposite limit (as far as we are investigating the properties 
of dense packing) . From a topological point of view, the densest circle packing is characterized by circles which are 
all in contact between neighbours. In such a packing the interstices between circles are all closed and surrounded 
by three circles. In this case, holds the identity 3F = 2E (any face has three edges and any edge divides two 
faces) and the vertex connectivity reach its maximum value: (z) Max = 6(1 — In general, less dense packing are 
characterized by lower value in the vertex connectivity (for example a typical connectivity for disordered packing 
of binary and polydisperse mixtures of circles is (z) ~ 3.75 [p7|). 



2.1 Apollonian covering 

An example of topologically densest circle packing is the Apollonian covering. In such a packing, any circle is in 
contact with 3 surrounding circles. A well known formula (The kiss precise | po| , gives a relation between the 
bends t\. = 1/Ri of the four "kissing" circles of radii R{, i = 1, ...4 

2{ef + e 2 2 + e 2 3 + e\) = (ej + e 2 + e 3 + e 4 ) 2 ■ (4) 



Following Coxeter 29 , one can consider the case in which the bends belongs to a geometric progression e„ = x v €q, 
where x is the ratio of the progression. Substituting into (||) one gets a sixth order equation in x with only one 
real root bigger than one: x ~ 2.89. In general, similar values for the ratio e^+i/e^ are characteristic of several 
covering sequences. For example, a value slightly less than 2.9, is the ratio at which converge the loxodromic 
spiral sequence J29|. These values corresponds to particular sequences selected from the whole covering set. The 
analytical evaluation of an average value of x corespondent to the complete set of covering sequences is still an 
open problem |32| . 

The Apollonian covering is a well-defined problem and constitute an excellent example of circle-covering. On 
the other hand, the aim of this work is the study of natural systems which have an higher degree of randomness 
in the packing strategy. 
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2.2 Disordered covering 

Let us consider a covering procedure that utilizes non overlapping circles with sizes distributed between two radii 
R-Max an d Rmin- Such a procedure starts filling the available space with the circles of higher radius Rq = Rhiax- 
Once the interstitial space between these circles become too small to contain circles with maximum radius the 
procedure reduces the radius of the new circles to insert R\ = Ro/x. The filling continues iterating this process 
gradually reducing the radii R v = x~ v Ro in order to insert circles of biggest radii compatibly with the interstitial 
sizes. Placing the centres in appropriate positions, this procedure generates the Apollonian sequences of tangent 
circles discussed in the previous paragraph. More generally, one can study a disordered covering where the circles 
are placed with the centres at random. Differently from the Apollonian case, in this covering the circles are not all 
tangent. In the Apollonian covering the interstitial space is always generated by three circles all in contact with 
each other. In the disordered case, the same three circles are -in general- not all tangent and the interstitial space 
is bigger. It follows that, in this case, the interstitial space can be filled with circles of bigger sizes than in the 
Apollonian case. Consequently, the ratio between the radii of the circles that form the boundary of an interstice 
and the circle that fills this interstice is lower in the disordered covering respect to the Apollonian case. Thus, the 
value x ~ 2.9 suggested as a limit of the ratio e^+i/e^ in the Apollonian sequences of tangents circles, represents, 
for the disordered case, an upper limit. 

It is now interesting to investigate the lower limit for x. This limit can be achieved by imposing the condition 
that the sum of the areas of all the circles introduced with the covering procedure converges to a finite value (which 
must be lower than the total area to cover). 

Consider, as before, a covering characterized by a sequence of radii in the geometrical progression e^ 1 = R v = 
x~ u R$. The number of interstitial circles introduced at any stage v is equal to the number of interstices of the 
system at that stage. This number is given by eq. @ 

*L=(^-l)K+X , (5) 

where V v is the number of circles at the stage v. 

The covering goes on, filling these interstices with F v new circles of radius R v +i. The number of interstices at 
the stage v + 1 is then 

F v+1 =(^-l S j{V v +F v )+ X ■ (6) 
Substituting the value of V v given by eq. (||) into (^) , one gets 

F v+l =(^f\F v , (7) 



that yields 

The total area covered by the circles at the stage v is given by the sum over the areas covered at any stage until 




A„=ir\V Ri + ^F i Rij=TrRi\V +F 0Tr ^y [- r ^)} . i<)i 

When v — > oo, this sum converge to a finite number only if (z)/(2x 2 ) < 1. This condition gives the lower limit 
for the sequence ratio 

*>Vt" • (10) 

Note that eq.(|^) have been obtained supposing the coordination number (z) and the ratio x being constant at 
any covering stage. This is not in general true. (The value of (z) is fixed and equal to 6 in the Apollonian covering 
of an infinite Euclidean plane but not in general. Moreover -also in this particular case- the ratio x is constant 
only in the asymptotic limit.) Should be noted that, the convergence condition for the total area, constrains only 



the asymptotic values of x and (z). (It is straightforward to see that, in the Apollonian case, equation (10) is 
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satisfied by the asymptotic values (x ~ 2.9 and (z) = 6).) In the asymptotic limit, the covering procedure is in 
general scale- invariant, as the only important metric parameter is the ratio between the sizes of the circles that 
generate an interstice and the size of the circle that fills this interstice. In this asymptotic limit the parameters (z) 
and x are related only to the covering strategy and are independent of the specific stage v. The inequality Jl(i| ) 
constrains these asymptotic parameters. 



3 Circle sizes distribution 

The circle-size distribution (N(R)) is related to the value of x, and to the number of cells introduced at any 
covering step. By associating the relation R u — R x~ u with eq.(||), one gets 

N(R„) = N , (11) 



where the exponent is 



, ■ (12) 

In a; 

Above we pointed out that the value of (z) is associated with the packing strategy: it is equal to 6 in the 
Apollonian case, whereas in the disordered covering it tends gradually to asymptotic values lower or equal than 
6. In the previous paragraph we discussed the bounds on the value of x. Substituting the lower bound given by 



eq.(10) one gets the upper limit for the exponent a 

b/<*> 
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a < V n -L = 2 = OCMax , (13) 



which is independent on (z). 

We discussed above that the upper bound of x is associated with the Apollonian covering. In this case we have 
x ~ 2.9 and (z) = 6 which leads to a = 1.03. 

The Haudorff Basicovitch fractal dimension of the packing can be directly derived from the size distribution 
one finds df = a. The analytical determination of df is a challenging problem of surprising difficulty. Only 
bounds are known 1.300197 < df < 1.314534 [|[ |j. The value here derived is quite far from the two exact bounds 
and from the value df = 1.305684 obtained by numerical calculations. It should be noted that our relation (JlJ) is 
an approximate solution and that the value of x = 2.9 that we assumed in ( [l3|) is only an estimation valid for a 
particular subset of all the possible covering sequences. On the other hand, values of 1 < a < 2 are in agreement 
with two exactly solvable models discussed in appendix [§] jl9| . In these particular cases (where x is constant) one 
have df = a = 1 and df = a = 1.585. 

Note that, in literature (for example |@, 0)) l ne s i ze distributions are discussed in the continuos limit. Here 
N(R V ) is a discrete distribution stating the number of circles with radius i?„. In the continuum limit, N(R V ) 
corresponds to the integral of the continuous distribution n(R) evaluated between -R^+i and R v . From eq.jll]) 
and using the identity R u +i = R v /%, one obtains n(R) — CR~ a ~ x with C — cxNoRq /(x a — 1) (n(R)dR being the 
number of circles with radius between R and R + dR) . 

3.1 Generalization to arbitrary dimensions 

It is possible to extend the results obtained in the previous paragraph spaces of any dimension. The calculus can 
be done following the same principal points utilized in the two dimensional case, extending the notions developed 
for the circles to D-spheres. 

Consider the network obtained by connecting with edges the centres of all the D-spheres mutually in contact. 
The number of vertices (V) of such a network is equal to the number of D-spheres. The number of cells (in the 
following indicated with C) is given by the relation 

C=^V , (14) 
no,D 
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where no : o is the mean number of D-cells incident on a vertex and no,D is the mean number of vertices bounding 
a D-cell. In two dimension, the incidence number 712,0 is the vertex connectivity (z) and 710,2 is the mean number 
of edges per face < n >. In a planar network these two parameters are connected trough the relation: 

(„) = 2{z) . (15) 

Substituting this relation into ( |l4| ) one gets immediately eq.(Q). For D > 2, 7t/j,o an d no,D are, in general, 
independent, free parameters related to the network peculiarities. ^ 

Consider the system at the covering stage v. At this stage, V„ D-spheres are packed. The number of interstices 
between them is C„ = (riD,o/no ; n)K- At the following stage u+1, these interstices are filled with C v D-spheres. 
Thus the total number of .D-spheres in the system become 

Vu+i = V v + C v = ( 1 + ^) V v . (16) 



Consequently the number of interstices is 



that gives 



n .D 



C v+l = [1 + 2M)C V , (17) 
no.D 



C u =[l + ^) C . (18) 
n ,D 



This equation generalizes eq.(|[). 

The covering procedure should be scaling invariant: in the covering strategy it is only important the ratio between 
the radii of the D-spheres that bound the interstitial space and the radius of the sphere that fills this interstice. 
The only important metric parameter that characterizes the covering procedure is the ratio x = R v /R v+ i. A 
limiting value for x can be obtained -as in the 2D case- imposing that the sum over the volumes, introduced with 
the covering procedure, converges to a finite value: 
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lim VftRf = Cbitf lim V \ = A . (19) 

This condition gives a lower limit for x: 



«D,0 ' 



x > l + ^ . (20) 

V n 0,D/ 

An upper limit for x can be calculated solving the algebraic equation obtained substituting e v with e v = x v ', 
into the generalization to any D of eq.(|j) |n|. (For example one obtain x ~ 1.88 for D — 3, x ~ 1.55 for D = 4, 
x ~ 1.39 for D — 5, x ~ 1.30 for D = 6 and a; = 1 for D — > 00). 



The size distribution can be obtained by associating eq.(18) with the condition R v+ i/R v = x. One has 



where the coefficient is 



ln(l + 

. (22) 



In 



2 For example, a particular case is the topological densest packing where any D-sphere is in contact with all the surrounding spheres. 
In this case any interstitial space is bounded by D + 1 D— spheres. In such a packing the network is made only by hyper-tetrahedron 
and one has tiq u = D + 1 (any hyper-tetrahedron has D + l vertices). On the other hand, the incidence number no,o is -also in this 
particular case— a free p ara meter which depends on the network characteristic. For D = 3, the maximum value of n^fi for the packing 
of equal spheres is 22.6 |J30|. By filling the interstices with lower sized spheres this number is reduced up to the lower limit of 12. 
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Figure 2: Circle-size distributions of two computer generated circle coverings, (a) Normal scale, (b) Double 
logarithmic scale. The full line is the average best fit. 



By substituting the lower limit of x (maximum packing, eq.(pfjj) ) one obtains the upper bound for the exponent 
ctMax = D. The lower bound for a can be calculated substituting into eq.(p2|) the maximum value of x (that 
corresponds to the Apollonian case) and the minimum value of the coefficient riD,o/ n o.D- This last parameter 
depends on the packing strategy. For D = 3, the Apollonian packing has: n^jj = 4, no.o > 12 and x ~ 1.88. In 
this case, one gets the limit of a ~ 2.2. 

The two models, discussed in app.[B| can be easily extended to the 3D case, one obtains a = 1.26 and a = 2 for 
the hexagonal and triangular models respectively. 



4 Computer simulations 

The covering procedure starts positioning the circles with maximum radius Rq = Ruax- The position of the centre 
of each circle is chosen at random. A new circle is added to the system only if it is non overlapping with any 
other pre-existing circle, otherwise another random position is chosen. This first step ends when all the interstitial 
spaces between circles have sizes smaller than the circle diameter and thus new circles with radius Rmux cannot 
be positioned any more. At the second step, the radius of the new circles is reducedi?i = Rq/x, the centre are 
chosen at random and a new circle is placed only if it is non overlapping with any other pre-existing circle. The 
procedure continues gradually reducing the sizes of the circles to place once the interstitial spaces become too 
small to contain new circles. No shifting or rearrangement are performed once a circle have been positioned. The 
procedure ends, after a finite number of steps, at the lowest radius R m in- 

Figures (^a,b) show the circle sizes distribution in normal and double logarithmic scale, for two simulations. 
The data have been obtained covering a rectangle of 600 x 450 pixels with periodic boundary conditions. The 
maximum radius was Ruax — 45 pixels and the minimum was R m in = 1.5 pixels. In accordance with the covering 
procedure, the circle size have been decreased between Rmux to R m in in 50 finite steps. In the simulations reported 
in fig.(|]), the total number of circles was respectively equal to 995 (squares) and 987 (triangles). 

The circle size distribution is characterized by one peak at the maximum radius and by a fast increment towards 
the minimum radius. The linear trend of the distribution, in the region of small radii (i? < i?Maa;/5), shown in the 
double logarithmic scale (fig.(Qb) ), suggests a power law N(R) cx R~ a . The best-fit estimation for the exponent 
give a = 1.84 ± 0.16 (squares) and a = 1.92 ± 0.19 (triangles), with confidence factors respectively equal to 93% 
and 92%. 

A power law for the sizes distribution was analytically predicted in the previous paragraphs and in appjA]. The 
best-fit values of the exponents a are in good agreement with the theoretical predictions for the case of disordered 
maximum circle-packing. 
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5 Structures formed by condensed tin droplets 



Figure (Q) is a SEM micrograph of Sn drops condensed on a hot, fiat alumina substrate. The system have been 
prepared evaporating the tin in high vacuum on a substrate heated at a temperature higher than the tin melting 
point p3| , p5[ | . After cooling at room temperature, one has a stable system of packed tin drops. The size-distribution 
of such a system was studied using SEM micrographs of magnifications 10 and 20 K ex. The micrographs were 
digitalized with a scanner. On each image an internal rectangular perimeter was defined and the drops inside and 
on the bound perimeter were encircled with circles of the same size of the surface occupied by the drop. The size 
distribution of these circles was studied. The boundary effects were taken into account counting 1/2 the circles 
crossing an edge of the perimeter and counting 1 /4 the circles containing a vertex of the bound-perimeter. In fig. (|^) 
the sizes-distributions for three different samples are reported in normal and double logarithmic scales. The three 
samples were deposited respectively with 2.5 g/m 2 (triangles), 2.8 g/m 2 (squares) and 3.2 g/m 2 (diamonds). The 
biggest drops have sizes of about lOOOnm, whereas the smallest (measurable) drop-sizes are 50nm. The number 
of drops encircled per each micrograph was about 2000-2500. The distribution is characterized by rather uniform 
drop sizes in the region of large radii and by a fast increment, as a power law, towards the smallest radii. 

Typically, one drop nucleates on a preferential centre and grows supported by condensation. This independent 
growth stops when the drop touch another surrounding drop. At this point, the two drops melt together and the 
formed new drop -eventually- melts with other surrounding drops in a chain of coalescence events. The coalescence 
leads to a new drop with volume equal to the sum of the volumes involved in the coalescence phenomena. This 
new drop occupies a surface lower than the sum of the surfaces occupied by the original drops. It follows that the 
coalescence mechanism liberates space on the substrate. On this surface the nucleation, growth and coalescence 
phenomena can start again. The final structure is therefore a consequence of a very complicated process which 
involve many body phenomena JlO| . Nevertheless the formation of the final structure can be investigated adopting 
a simple point of view. 

The system of quenched drops showed in fig.(^|) is apparently similar to the systems studied, analytically 
and numerically, in the previous paragraphs. This similarity is not only apparent: the nucleation, growth and 
coalescence mechanism leads to systems where the drops occupies the maximum possible surface and have the 
biggest sizes, compatibly with the coalescence that prevents the drops overlapping. The system evolves through 
self-similar configurations where the size distribution is, at any time, the maximum packing distribution. 

The similarity between this two systems is confirmed by the drop size distribution: in the region of little radii 
the distribution follows a power law N(R) cx R~ a with exponent a ~ 2. The best-fit give a = 1.96±0.13 (triangles 
in fig.(D) ), a — 2.05 ± 0.11 (squares) and a — 2.1 ± 0.23 (diamonds), with confidence factors respectively equal to 
94%, 95% and 87%. These values are in accordance with the analytical predictions and strongly suggests that the 
system morphogenesis is ruled by the mechanism of the maximum circle-packing. 

In the region of big radii (R ~ Ruax > 500 nm) the drops have a rather uniform sizes. This region of uniform 
sizes is the memory of the first stage of the deposition where the drops grow independently. Indeed, in this stage the 
relative volumes of the drops are proportional to the sizes of the Voronoi cell constructed around the drop centers 
and a rather uniform distribution is expected. In fig.(||) this uniformity at R ~ Rmux is not particularly 
manifest since the number of drops in this region is little in comparison with the number of small drops ( for 
example in a sample deposited with about 3 g/m 2 , one has about 10 drops per [im 2 of about O.l^tm and 1 drop 
per /im 2 of l^m). Despite the fact that their number is little, the big drops cover a large amount of surface and 
contain the main part of the volume deposited (using the same example reported above, it is straightforward to 
see that the big drops cover an area 10 times larger and occupy a volume 100 times larger than the small drops). 
Note that, a peak at Ruax is present also in the computer simulations. A deviation from the maximum packing 
law in the region close to Ruax is not in contradiction with the theoretical predictions. The analytical result 
N(R) cx 1/R a , concerns only the asymptotic limit of the distribution [y — > oo i.e. R/Rhiax — ► 0) where the system 
is supposed to be scale-invariant. 

Studies reported in literature |nj [ll], |l2|, [13], [u], [l5| (L6, 17 , show that the system of condensed drops evolve in 
time increasing the drop sizes, decreasing the drops number and maintaining constant the coverage. In the systems 
of tin drops here studied, the coverage have been experimentally found in the range 60% to 65% |2^ ], p5| . A 
small value of the coverage (typically in breath figures is equal about 55% |l5|]) is a consequence of coalescence that 
continuously liberates space. Consider a configuration of three external drops of radius Rq and one interstitial drop 
of radius R\. The collapse of these four drops into one trough coalescence liberates a certain amount of surface. 
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Figure 3: SEM micrograph of a "breath figure" generated by tin drops deposited on a hot alumina substrate. 
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Figure 4: Size distribution of the tin drop diameters for three different samples, (a) normal scale; (b) double 
logarithmic scale. In ordinates, N(R) is the number of drops per fxm 2 . The full line is the average best fit. 
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The fraction of area liberated (Area occupied by the drops after coalescence/ Area occupied before coalescence) is 
slightly dependent on the packing strategy. It is equal to 32% if the four drops are all "kissing" each other and 
equal to 35% in the case of disordered maximum packing (Rq/Ri = x ~ v3). The presence of other interstitial 
drops that participate to the coalescence phenomena reduces the previous values within 1%. Starting with an 
initial -hypothetical- coverage of 100%, the coalescence reduces the coverage to 68% and 65% in the two cases 
discussed above. These values are close to the coverage experimentally observed. Note that these values are scale 
invariant and depends only on the local configurations of drops (i.e. on the packing strategy). The scale invariance 
implies that the fraction of area liberated by coalescence is independent of the drops radii (i.e. independent of the 
amount of tin deposited) thus constant during the deposition. 



6 Conclusions 

Circle-covering is a good model-system for many two dimensional natural cellular system where the plane is filled 
in the most efficient way compatibly with the cell shapes and sizes. It has been found (jQ §|| and app.[A|) that 
polydisperse circles packed in a dense way, have a size distribution that follows the general law N(R) cx R~ a . 
This distribution is a consequence of the scale-invariance in the packing strategy. The range of variability of the 
exponent have been calculated. We have obtained a maximum value of a.Max = 2 for the general disordered case, 
where the circles are arranged at random following the only constraint of non overlapping. The minimum value 
has been estimated a m i n ~ 1 and corresponds to an approximate solution for the Apollonian packing (§ [2.l| and 
§||), and to an exact solution for the "hexagonal" filling model (app.[B|). 

These results, obtained for two dimensional circle packing, have been extended to packing of spheres in spaces 



of arbitrary dimension (§3.1). We found that, in D-dimensional spaces, the size distribution of densely packed 
D-spheres follows -as in 2D- a power law N(R) cx R~ a . In this generalized case, the maximum value of the 
exponent is ot-Max — D and is associated with disordered packings. The minimum value is associated with the 
packings of tangents D-spheres. For D = 3, this value have been estimated equal to a = 2.2 for the Apollonian 
case, whereas for the hexagonal filling model we found a = 1.26. 

Computer simulations of two-dimensional circle-packing confirm the analytical predictions: in the disordered 
case the size-distribution follows a power law with exponent a ~ 2 (see §0) . 

The formation of breath figures have been interpreted in terms of circles packing regulated by coalescence (§||). 
The mechanism of formation of this figures is very complicated since it is a dynamic system where the evolution 
involves many-body phenomena. On the other hand, the structures formed by the drops instant by instant have 
strong similarities with the structures generated by packing circles. In particular, the size distribution in the 
region of little radii, follows the same power law as the disordered, dense circle-packing with exponent a ~ 2. The 
coverage, evaluated through the fraction of area liberated by coalescence in a system of drops that follows the 
maximum packing distribution, is in agreement with the experimental observation. 
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A Size-distribution and covering strategies 



In the second and third paragraphs we found that the size distribution of densely packed circles follows a power law 
(cq. ( pl|) ) where the exponent a is related to the topological properties of the Dodds network and to the parameter 
x (eq.(|12[) ). This is a general result where the only a-priori hypothesis is the existence of a Dodds network with 
convex cells (this condition can be considered as the definition of dense packing). In this appendix we generalize 
these results using an approach which doesn't needs the definition of the Dodds network. 

Consider a covering procedure which fills the available space with non overlapping circles starting from the 
circles of bigger sizes Ruax and then gradually reducing the sizes R v < Ruax in order to fill the interstitial spaces. 
The number of circles introduced at any covering step and the relative sizes depend on the covering procedure. In 
the second paragraph we pointed out that the scaling condition implies that the number of circles inserted at any 
covering step grows as a geometrical progression N v cx a v and the size of the circles decreases as a geometrical 
progression R v oc Rq/x" . As a consequence the size distribution results a power law N(R U ) cx R~ a with exponent 

Ina . . 

a = — . 23 
mx 

In the case of dense packing (i.e. when it is possible to define a Dodds network with convex cells) we found 
a = (z)/2 (eq.@ ) and the bound x > y/(z)/2 (eq.© ). 

More generally one can construct the Delaunay triangulation with the centre of the circles as vertices. This 
triangulation is always well defined and (for an infinite system) the average connectivity is (z) = 6. Different packing 
strategies differentiates in the number, the sizes and the positions of the new vertices (i.e. the new circles) to insert. 
From a topological point of view, there are only two position where the centre of the new circles can be placed: 
inside a triangle or on an edge (the centre on a vertex is forbidden by the condition of non overlapping). Suppose 
that the packing procedure places a new circle with probability p inside a given triangle and with probability q on 
a given edge of the Delaunay triangulation. If V v is the number of vertices (i.e. of circles) at the covering stage v, 
it is easy to prove (following the same arguments of §g) that at the next stage this number is 

V v +i = (1 + 2p + 3q)V v . (24) 

From relation (p4|) follows immediately that the number of circles grows as a geometrical progression (N v = V v oc 
a v ) with coefficient a — (1 + 2p + 3q). 

The argument used in ^ to find the lower bound on the parameter x (eq. ( [To| ) ) can be directly extended to the 
present case. One gets x > ^/a, which, substituted in (p3|), gives the upper limit an ax = 2. 

The Apollonian covering is a particular example of the covering procedure here discussed. In this case one has 
p = 1, q = and x ~ 2.9, which gives a ~ 1. 



B Hexagonal and triangular filling models 

Following the work of Bidaux et al. |l9[ ] , let us describe two simple circle-covering problems where the exponent 
a can be calculated exactly. Originally these models have been proposed as simplified geometrical problems to 
estimates two bounds for the fractal dimension in the Apollonian packing. In the present work these two models 
are presented as good examples of circle-covering which do not necessarily corresponds to approximations of the 
Apollonian case. 

Let us start with one triangle. In the first model (triangular), one insert a new triangle with vertices in the 
mean point of the edges of the original triangle. In this way the original triangle is divided in 4 identical triangles 
which arc similar to the original. Then, one insert a circle inside the central triangle and iterate the procedure on 
the three external triangles. At the beginning one starts with 1 circle, at the first step one insert 3 new circles, at 
the v th step the number of new circles inserted is 3". The parameter a (see. app.|X|) is then equal to 3. The ratio 
x can be easily derived by observing that the triangles are all similar and that at each step of the sequence the 
sizes of the edges are reduced by a factor 2. Consequently one has x — 2. Substituting the parameters a and x 
into eq.(p3|) one gets the exponent a — In3/ln2 ~ 1.58. 

In the second model (hexagonal) one inscribe first an hexagon with vertices which are dividing in three equal 
part the edges of the original triangle. Then a circle is inscribed inside the hexagon and the procedure is iterated 
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on the three remaining external triangles. As before one obtains a — 3, whereas the scale factor is x — 3. By 
substituting into eq. (|23|) we get the exponent a = 1. 

Note that, if one use equilateral triangles the coverage 9 (Area covered by the circles/ Area of the original 
triangle) can be calculated exactly. One has = 7r/(3V3) — 0.6 for the triangular model and 9 — 7r/(2-\/3) ~ 0.9 
for the hexagonal. These value can be increased by inserting new triangles in the free interstices between the circles 
and iterating the procedure. 
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